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ABSTRACT 



The dynamical evolution of super star clusters (SSCs) moving in the back- 



\ ground of dark matter halo has been investigated as a possible event causing 



the observed correlation between the mass of galactic bulge, Mtuige, and the 
mass of its central black hole, Mbh- The involved physical processes are the 



o 

Qh' sinking of SSCs due to the dynamical friction, and the stripping of SSCs on 

O \ their way to the center. Model calculations show that only sinking of circum- 

C/3 . nuclear SSCs contribute to both the growth of the central object and the 

■ 

formation of the galactic bulge at the early stage. On the assumption of a 



universal density profile for the dark matter halo, and an isothermal model 



H ' for the SSCs, our simulations have yielded the mass ratio of the central objects 

to the bulges formed this way to be about a few times 10^'', less than the 
observed median value for early type galaxies. It is, however, consistent with 
the observed mass ratio for disk spirals, implying that the proposed scenario 
might be a possible event for the formation of bulges and central black holes 
of late type galaxies and for Mbh — Mbuige correlation of disk galaxies. 

Key words: black hole physics - Galaxy: center - Galaxy: bulge - Galaxy: 
kinematics and dynamics 



1 INTRODUCTION 

Recent observations with unprecedented high resolution have firmly established that many 
galaxies, whether active or not, host supermassive black holes (SMBHs) in their centers (Ko- 
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rmendy & Richstone 1995, Kormendy & Gebhardt 2001). These observations have revealed 
correlations of the mass of the SMBHs to the luminosity (or mass) of their bulge (Kormendy 
& Richstone 1995, Magorrian et al. 1998, Wandel 1999), or to the luminosity-weighted line- 
of-sight velocity dispersion within the effective radius with much smaller scatter (Gebhardt 
et al. 2000, Merritt & Ferrarese 2000). However, no correlation has been found between 
the BH masses and the total luminosities of host galaxies (Kormendy & Gebhardt 2001). 
The median black hole (BH) mass is 0.13% of the mass of the bulge (Kormendy & Geb- 
hardt 2001) for early type galaxies, or much smaller for late-type spirals (Salucci et al. 2000, 
Kormendy & Gebhardt 2001, Gebhardt et al 2001). 

The latest observations with X-ray Observatories provide evidence showing that the 
formation of black holes is connected with starbursts (Matsumoto et al. 2001, Fabbiano 
et al. 2001, Strickland et al. 2001). Matsumoto et al. (2001) and Strickland et al. (2001) 
argued that there exist intermediate-mass black holes (MBH) in the off-nuclear compact 
X-ray sources in M82 and NGC3628, respectively. The formation of SMBHs through these 
MBHs has been explored soon after (Ebisuzaki et al. 2001). Haenelt & Kauffmann (1999) 
have also studied the relation between SMBHs and the formation of galaxies. On the other 
hand, a beam model has been proposed (King et al. 2001) to ease the difficulties that the 
unbeamed models run into, where MBHs are required. 

The correlation between the mass of SMBHs, Mbh, and the mass of the bulge, Mbuige, 
implies that the growth of the central BH and the formation of the bulge are probably 
caused by a same physical process, or event (e.g. Kormendy & Gebhardt 2001). This "event", 
however, has not been clarified yet. 

In this paper, we report our study on the dynamical evolution of circumnuclear SSCs 
moving in the dark matter halo as a possible "event" for the observed relation of Mbh to 
Mfjuige- We find that the sinking of the SSCs, along with the tidal stripping on their way to 
the center, could contribute to simultaneous growth of Mbh and Mbuige- Based on a set of 
reasonable parameters, the mass ratio Mbh /Mbuige derived from our numerical simulations 
coincides with the observed values of disk spirals. 

2 MODELS 
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2.1 SSC model 

Recent observations of starburst galaxies with high resolution have revealed that there are 
many compact, young and very luminous SSCs in the central regions of galaxies (Shaya et 
al. 1994, Surace et al. 1998, Surace & Sanders 1999, Whitmore et al. 1999, Scoville et al. 
2000, de Grijs et al. 2001). Numerical simulations have also shown that interaction or merger 
among galaxies can trigger strong starbursts around the nuclear regions. While galaxies are 
merging or interacting, some gas components could lose their angular momentum and fall 
into the central region (Mihos & Hernquist 1996, Barnes & Hernquist 1996). Then, rather 
high pressure of the warm interstellar gas could induce global collapse of giant molecular 
clouds and thus forming circumnuclear SSCs (Jog & Solomon 1992, Harris & Pudritz 1994). 

On the one hand, there exist some observational facts about SSCs, but not many con- 
clusive results about their dynamical properties are reached so far. On the other hand, the 
SSCs are believed to be the progenitors of present-day globular clusters (GC) (e.g. Smith & 
Galagher 1999, Origlia et al. 2001). It would then be reasonable to model SSC by using both 
clues mentioned above. We assume that SSCs have a similar mass spectrum as the initial 
GC mass function but with larger mean value. According to Vesperini (2000,2001), we take 
the following log- normal mass function for SSCs (SSCMF), with its mean at 5 x lO^M© 
following the investigations on SSCs (e.g. Surace & Sanders 1999, Origlia et al. 2001, Surace 
et al. 1998, Cen 2001, Meylan 2001), 

Logio{M) ~ N{exp{mean) = 5 x 10^, variance = 0.08) (1) 

With this SSCMF, we generate randomly 100 sets of SSCs, each containing 100 SSCs (for 
taking 100 SSCs, see Combes 2001, Fellhauer 2001, Whitmore et al. 1999). 

For simplicity, an SSC is modeled as a truncated isothermal sphere with three parameters: 
the central density (p^), the velocity dispersion (a) and the initial truncated radius (Rq). 
Taking a lower and an upper limits of pc to be 5.3 x 10^Mq/pc^ and 3.4 x from 
observed ones (Larsen et al. 2001, Campbell et al. 1992), we assume a linear function of 
cluster mass, M, for pc{M). And with Rq taken to be the local tidal radius, a can be derived 
from M and pc- The resulted a also increases with M. 

The initial distribution and evolution of SSCs in a starburst galaxy is not clear, though 
it is very important for the problem we are investigating. In our simulation, all of the 100 
circumnuclear SSCs in each set are initially placed at a distance of Ikpc from galactic center, 
typical locations for circumnuclear SSCs, and assigned with local circular speed. As will be 
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discussed in the last section, we also tried larger distance and found that more distant SSCs 
are generally irresponsible for the formation and early growth of bulge and central compact 
objects. 

2.2 Background 

As we are considering the formation process of bulge and central massive objects, the dark 
matter dominates not only globally, but also in the central region at early stage. Therefore, 
only dark matter halo is considered initially in our simulation. For the dark halo, we assume 
the universal density profile (Navarro, Frenk & White 1997), which can be written as (e.g., 
Binney et al., 1998) 

pf^y^) = -( — (2) 

where r is the distance from the halo center, and, an and M^h are two parameters. The value 
of ah is connected with the extension of the halo and it may change with time. However, how 
it changes is not well quantified. So, we take two rather different values of ah for comparison. 
Later, we will see that the results are not too much different, and so, acceptable. Besides, 
since we are concerned with the formation and the early growth phases of galactic bulge 
and its massive central objects, of which the time duration is taken to be IGyr, ah might be 
well approximated as time-independent. Following El-Zant et al. (2001), we take ah = 6kpc 
and lOkpc, and, Moh is derived from ah with the condition that there are lO^^M© interior 
to 200kpc. 

Our test calculations show that the first fallen SSCs, the massive ones in each set of SSCs, 
can contribute about 2 x IO^Mq to the central globe of radius Ipc, while in the same region 
there are dark matter of less than 10^ Mq. In the meantime, the mass in the innermost region 
(interior to a few lOpcs from the halo center) can be significantly changed by the fallen SSCs, 
although outside about 200pcs remains dark halo dominating. The background variation of 
this kind has substantial effect on the further mass contribution to the central region from 
stripped SSCs. In order to account for the effect of fallen SSC mass, which is characterized 
by a steeper mass density profile near the center than the universal one, a truncated singular 
isothermal sphere is added after the most massive SSC has fallen. 
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2.3 Dynamical friction 

In the following, M and Vm (with Vm = \Vm\) denote, respectively, the mass and velocity 
of cluster experiencing the dynamical friction. Assuming a Maxwellian velocity distribution 
with dispersion abackground of background matter, composed of particles with mass much 
smaller than M, the dynamical friction formula writes (e.g., Binney et al, 1987) 

^t7 27rlog(l + A2)(:;2^p 2X . ...^-^ . . 

dVMjdt = ,,3 ^[erf(A:) - — exp(-X')]VM (3) 



v^j ' ' ' 

where erf is the error function, and, 

h V2 



< 



GM 

(4) 



ackground 

The quantity h^ax is the so-called maximum impact parameter and Vtyp a kind of typical 
module of the relative velocity between M and a background particle. Neither hmax nor 
Vtyp is precisely defined. Fortunately, uncertainty in either quantity causes no significant 
difi^erence in the resulted values of the dynamical friction. Following Binney et al. (1987), we 
use hmax = ^kpc and take Vtyp = Vm- The velocity dispersion Chackgroundir) can be roughly 
estimated from the Jeans equation. 

2.4 Stripping 

It is assumed that the stellar mass outside a sphere, the radius of which is denoted as 
Rt, approximating instant Hill stable region around the SSC center, will be stripped. Since 
the stripping is processed continuously as r (the distance between halo and SSC centers) 
decreases, only a thin outer layer is to be stripped at a time. Therefore, in an average sense, 
the stars stripped when the SSC goes from r to r — dr are contributed to a region radially 
bounded by r + Rt{r) and r — dr — Rt{r — dr). As a first-order approximation, the mass of the 
stripped stars are considered to be, at some later epoch, uniformly distributed in the shell 
bounded by r+Rt{r) and r—dr—Rt{r—dr). By summing up all of the stellar mass stripped at 
various rs, the stripped stellar mass distribution can be derived. The SSC's mass contributed 
to the galactic center (taken as a globe with radius Ipc) is just the remaining mass of the 
SSC when r -|- Rt{r) < Ipc plus the previously stripped mass inside the above-mentioned 
Ipc globe. 

Obviously, if a massive single object is embedded in the center of SSC, stripping cannot 
be proceeded further when only this object is left. X-ray observations discovered a lot of the 
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Figure 1. Mbh vs Mi,.^ig^. Legend: solid line — results for early type galaxies; ★ — observational results for Sb-Im galaxies; 
open and filled circles — simulation results with and without non-striped compact cores (or BHs) embedded in the SSCs, 
respectively. 

so-called super-Eddington sources associated with SSCs (Matsumoto et al. 2001, Strikland 
et al. 2001). However, whether there are MBHs, ranges from several hundreds to about one 
thousand solar mass (Ebisuzaki et al. 2001), or the observations are only due to beam effect 
(King et al. 2001) is still not clarified. If massive black holes do form in SSCs, they would 
most likely be at the center of the SSCs. Therefore, in our simulations, we consider two 
extreme cases: the stripping is not allowed when the mass of the stripped SSC is less than 
IMq and lOOOM©, respectively. 

3 RESULTS 

The results about the relation between SSCs' mass contributions to Mi,uige and Mbh are 
summarized in Fig.l, where Mipc, Mioopc and M2oopc are the mass contributions of the SSCs 
to the globes of radii Ipc, lOOpc and 200pc, respectively. The observational results for early- 
type galaxies (Kormendy & Gebhardt 2001) and six Sb-Im galaxies (Salucci et al, 2000) are 
also shown. In our simulations, 94 - 96 percent of the SSC sets, each representing SSCs in 
a single galaxy, can contribute mass to the inner Ipc globe. Here we take two values, lOOpc 
and 200pc, for the radius of bulge. As shown in Fig.l, no substantial differences between 
these two cases. Noting further that the stripped SSC mass interior to lOOpc is generally no 
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more a small quantity in comparison with the dark halo mass in the same region, to take 
lOOpc as the radius of a bulge (or pseudo-bulge) form this way might be appropriate. 

By comparing the left and right panels of Fig.l, one finds that Mhuige — Mbh relation 
has no fundamental disparity for different values of ah- This implies that similar results 
could be obtained for time-dependent ah, regardless how it changes with time. Besides, the 
results for SSCs with and without central black holes are also similar to each other. This 
is because the mass of the formed central object of galaxy, much larger than that of the 
assumed SSCs central black hole {10^ Mq), is mainly comes from the most massive SSC As 
a result, the obtained Mbuige — Mbh relation might be valid for generic very-late type disk 
galaxies, provided that the bulge-BH do form in this way. 

As can be seen from Fig.l, Mbh /Mhuige in our simulation is smaller than that of the 
observed one for early type galaxies. And, as already stated, strong tidal force of galactic 
central object previously formed from the most massive SSC prohibits the non-single ob- 
jects sinking into the Ipc globe. This implies that the bulge mass will be increased faster 
afterwards, and so, even smaller value of Mbh / M^uige is expected at some later stage. Our 
results, however, consistent with the observational ones for disk spirals (Salucci et al. 2000, 
Kormendy & Gebhardt 2001, Gebhardt et al. 2001), implying that the proposed scenario 
of Mbh — Mhuige formation might be valid for very-late type disk galaxies. Possibly, the 
Mbh — Mbuige correlation for galaxies of various types of Hubble sequences might not be 
linear, which is a possibility discussed in the case of M33 (Gebhardt et al. 2001). To detect 
the BH masses of less than IO^Mq in galaxies would be crucial for deliberating the linearity 
of the Mbh — M^^ige correlation. 

4 DISCUSSION 

In our simulation, we have made some assumptions either due to lack of knowledge or for 
simplifying the simulations. 

An important assumption is that the background is composed of dark matter only, though 
we have considered the background variation later. Indeed, this assumption is compatible 
with what we investigate in this paper — the formation of bulges and the growth of central 
black holes at the early stage. In accordance with this, the SSCs we have studied are cir- 
cumnuclear ones assumably originated from the mergers of very late-type galaxies. These 
galaxies are sources with disk components only. In this case the dark matter of a few times 
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IO^Mq dominates over luminous systems inside 1 kpc from the center. As a first step of our 
investigations, it would be reasonable to make such an assumption. On the other hand, it 
wiU be interesting to see what the theoretical correlation will be if more components are as- 
sumed, e.g. disks and a pseudo-bulge, in the background for our simulations. In other words, 
the logical, next step of our investigations is to study a way for the first formed bulge and 
the central BH to grow further, and to see what the ratio of Mbh / ^huige "will be. Probably, 
a new merger is needed. That is, a new merger occurs between two disk galaxies with small 
central BHs and pseudo-bulges. A study of this kind is under our consideration. 

Another point we have not considered in our simulations is the distance distribution of 
SSCs, on an average of Ikpc for circumnuclear SSCs. SSCs at different distance would make 
different contribution to both bulge and central object masses. M33 might be an example 
in favor of considering distance distribution for SSCs. A compact star cluster is observed 
to locate at its nucleus (e.g. Gebhardt et al. 2001), which is more like a GC based on its 
dynamical parameters (Kormendy & McClure 1993). The formation of this GC might find 
its way in the scenario proposed in this paper. 

Besides, in our simplified treatments, we include no effects of non-spherical SSCs mass 
distribution, of background rotation, and of others. We'll cooperate these effects in our future 
studies. 

The survival of SSCs over Gyrs is a key question for the scenario presented in this 
paper, which is also crucial for the hypothesis where the SSCs are progenitors of present- 
day GCs. This question has been analyzed by Origlia et al.(2001) for one SSC in NGC 
1569, NGC1569-A1. They found that this SSC has a standard Salpeter initial mass function 
with no truncation at lower mass limit, implying that it can evolve into a system similar 
to the present-day GCs. Forbes, et al.(2001) have claimed recently that the formation of 
a bulge/spheroidal stellar system is accompanied by the formation of metal-rich GCs (red 
GCs). Adopting the scenario of SSCs as the progenitors of GCs, our proposed processes for 
the formation of bulge, i.e. through sinking of SSCs, along with the tidal stripping, thereafter 
might be a way for their claim. 

We have also performed numerical simulations for the SSCs located far away from the 
center, without considering the effect of tidal stripping. The results show that these SSCs 
with masses of 10^ — 10^ Mq will stay in the external regions. Cen (2001) proposed that 
the external SSCs are formed from gas-rich sub-galactic halos triggered by the reionization 
of the universe. These young stellar systems with masses of 10'^ — IO^Mq are suggested 
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as progenitors of the present-day halo GCs (blue GCs). Our study on the external SSCs 
supports Gen's suggestion in the sense that they will stay well outside the galactic central 
regions. 
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